LTH 580 



Three loop anomalous dimension of the second moment of the 
transversity operator in the MS and RY schemes 

J. A. Gracey, 
Theoretical Physics Division, 
Department of Mathematical Sciences, 
University of Liverpool, 
P.O. Box 147, 
Liverpool, 
L69 3BX, 
United Kingdom. 



Abstract. We compute the anomalous dimension of the second moment of the transversity 
operator, ■\\)a iiv D p tp, at three loops in both the MS and RI' schemes. As a check on the result we 
also determine the 0(1/Nf) critical exponent of the n-th moment of the transversity operator 
in d-dimensions in the large Nf expansion and determine leading order information on the n 
dependence of the anomalous dimension at three and four loops in MS. In addition the RI' 
anomalous dimension of the non-singlet twist-2 operator, ^^D u ip, is also determined. 
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1 Introduction. 



Recently there has been renewed interest in the transversity distribution of the partons in the 
nucleons. Introduced originally in ^ [21 El the transversity measures the difference in probabilities 
of finding a quark in a transversely polarized nucleon which is polarized parallel to that of the 
nucleon with that in the antiparallel polarization. Although there has not been as large an 
experimental activity in this area compared with the conventional distribution functions of deep 
inelastic scattering, it is possible RHIC may take transversity data in the near future, [I]. (For 
recent reviews see, for instance, pHE]-) Therefore, to improve our interpretation of the results it 
is necessary to extend our theoretical understanding to a similar level as that achieved in deep 
inelastic scattering. For instance, there has been a large activity in determining the three loop 
anomalous dimensions of the twist- 2 unpolarized and polarized flavour non-singlet and singlet 
operators as a function of the operator moment, n, IH1 E3 El EEH • These three loop results 
are necessary for the full two loop renormalization group evolution of the structure functions. 
Exact analytic results as a function of n are required since the Mellin transform with respect 
to n determine the parton splitting functions as a function of the momentum fraction, x. Such 
three loop calculations, however, require a huge increase in resources compared with the earlier 
one and two loop results of ^21 El El El El El El EES]- For instance, due to the presence 
of a large number of Feynman diagrams with and without subgraphs, one encounters nested 
or harmonic sums whose mathematics has had to be studied and developed in, for example, 
|2U| . Moreover, this has then to be encoded in the symbolic manipulation language Form, 
|21j . in order to handle the incredibly large amounts of algebra. Whilst results for non-singlet 
operators are expected soon it is clear that the same machinery can be applied to determining 
the anomalous dimensions of the moments of the transversity operator as they differ only in their 
7-algebra structure. For instance, the transversity operator involves the antisymmetric tensor 
a nv _ lpyi^Kj i ns t eac l f one ^-matrix. Whilst an analytic result for transversity at three 
loops is some way off, one can turn to the original three loop approach of [7j and determine 
information on the low moments of the anomalous dimension of the operator. Indeed in [22] the 
anomalous dimension of the first moment is available at three loops in QCD, building on the two 
loop result for this tensor current, [2H]. Moreover, the anomalous dimension of ^a^ u tp is also 
available in QED at four loops in MS in the quenched approximation, |2l]. In this article we 
will extend [22] by computing the anomalous dimension of ^pa^DPip at three loops where is 
the usual covariant derivative. Although one motivation for the result lies in the provision of an 
independent calculation which will be useful to check against the result as a function of n when it 
is determined, there are several other reasons for concentrating on this one operator. These lie in 
other theoretical techniques to determine information on the associated matrix element. Using 
a lattice regularization one can deduce numerical estimates for matrix elements. However, such 
results are necessarily in a lattice regularization scheme and need to be matched to the standard 
MS renormalization scheme. Therefore, we will not only determine the anomalous dimension at 
three loops in MS but also in the RI' scheme which denotes the modified regularization invariant 
scheme, |25j . This scheme was introduced in |25j and developed for the problem of improving 
estimates for quark masses in |26| l27j. For applications to similar problems in applying the lattice 
to construct matrix elements relevant in deep inelastic scattering see, for instance, |281l29j . It is 
also worth drawing attention to an alternative lattice approach to constructing moments of the 
parton distribution functions, |3U| . This uses the Schrodinger functional technique of |31| and is 
non-perturbative. More recently the RI' scheme has been examined in detail in jH2j where QCD 
has first been renormalized at three loops before reproducing the results of |27| at three loops for 
the scalar current. Hence the RI' expression for the anomalous dimension of the tensor current 
was determined to the same order, [221 • By contrast to MS the anomalous dimensions of these 
currents cease to be independent of the covariant gauge parameter. However, only results are 
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required in the Landau gauge for the lattice matching. Nevertheless from the point of view of 
internal consistency in multiloop computations we will perform our calculations in an arbitrary 
covariant gauge before deducing the requisite conversion functions in the Landau gauge. The 
disadvantage of an arbitrary gauge calculation is that the computation is algebraically more 
intense which is further complicated by the fact that the presence of a covariant derivative, 
itself, introduces another level of integration by parts, also slowing the procedure. Therefore, 
we have chosen to only focus on the second moment. Choosing the Feynman gauge would speed 
the computer algebra manipulations but would not be useful for the necessary matching which 
is gauge dependent. Further, matrix elements for higher moment operators from the lattice are 
currently much harder to extract. 

To be more specific about the results we will report here, we will renormalize both the 
(flavour non-singlet) operators tp^^D v ip and ip(j^ v D p ip at three loops in MS and RI'. Although 
the anomalous dimension of the former operator is known in MS, |12[ llH [HH ll5 [nT)} ll7 [llM[ ll9 [ l7]. 
we will use it as a check on the computer algebra programmes we have written in the symbolic 
manipulation language Form, [2J|) before extracting the RI' anomalous dimension. Then the 
Feynman rule for the operator insertion in a quark two-point function will be replaced with that 
for the transversity operator. One check on both the MS results will be that the anomalous 
dimensions of both operators will be independent of the covariant gauge parameter, a. The 
method of calculation is to apply the Mincer algorithm, written in Form, [HI]) which 
determines the poles in e in dimensional regularization, with d = 4 — 2e, as well as the finite 
part of massless three loop two-point functions. Each operator will be inserted at zero external 
momentum to form a two-point function. Moreover, lattice calculations are for the operator 
inserted at zero momentum so the renormalization scheme conversion functions are required for 
this configuration. We have used Qgraf, 35 , to generate the Feynman diagrams to three loops 
and the output has been converted to the format used by the Form version of Mincer. 

The paper is organised as follows. In section 2 we introduce the method we will use by 
considering the twist- 2 flavour non-singlet operator and reproduce the three loop anomalous 
dimension for the second moment in MS before deriving the RI' scheme result. This approach 
is extended in section 3 to the case of the second moment of the transversity operator. The 
conversion functions which allow one to translate from one renormalization scheme to another 
for each of the operators we consider are discussed in section 4 and provide a partial check on 
our results. As an independent check on our three loop result for the transversity we compute 
the large Nf leading order critical exponent for the nth moment of the operator using the large 
Nf critical point technique in section 5. Finally, we conclude with discussion in section 6. 



2 Second moment of non-singlet twist-2 operator. 

As we are considering operator insertions similar to the tensor current which was studied in 
|22l I32j . we briefly recall the RI' renormalization scheme definitions. To determine the correct 
finite part of the Green's function with the operator inserted we must project out all possible 
components consistent with Lorentz symmetry. For instance, for the flavour non-singlet operator 
ip^D^ip we have in momentum space, 

= 4%>)(7V + tV-^t) 

+ Ef m {p)Up^-^A (2.1) 
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where p is the momentum flowing through the external quark legs and S denotes both the 
symmetrization with respect to \x and v and that the operator is traceless*. We could have 
introduced the null vector, A„, to achieve this but it is not necessary here as we will in fact 
reproduce the known anomalous dimension for tp^/^D v ip. To be more specifc 



a 



(2.2) 



The various components are determined through 



S (2) (n) 



1 



8(d- 1) 



(p) 



4(d-l) 



2tr l^ji--jjv 1 G^ illD ^(p) 



(d + 2) tr I I p^ - ) Gj? (p) 



d 



(2.3) 



To renormalize (^(p) [tv^^-D^^O) 4>(~P)) m our symbolic manipulation approach we follow 
the method of jSHl- There the bare Green's function is computed in terms of the bare coupling 
constant and covariant gauge parameter. Renormalized variables are introduced through the 
usual renormalization constant definitions such as g = n e Z g g where g is the coupling constant 
appearing in the covariant derivative = + igT a A^, \x is the renormalization scale intro- 
duced to ensure the coupling constant remains dimensionless in d-dimensions and T a are the 
generators of the colour group whose structure functions are f abc . The renormalization constant 
associated with the operator, Zq, is deduced by ensuring that the Green's function is finite after 
multiplying by ZqZ^ where the quark wave function renormalization constant, Z^, arises from 

the external quark legs. For the RI scheme it is the first component which determines Zq 
since this renormalization constant is defined by, |251 132j . 



lim 



7 Rl' 7 Rl' y(l) ( s 
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(2.4) 



In other words, unlike the MS scheme in addition to the poles in e the 0(1) piece is absorbed 
into the renormalization constant. Though it is important to note that this renormalization 
constant will also render the other components finite but not necessarily zero or unity. 

Having summarized our method of calculation, we now record our results for the renormal- 
ization of tp^^D v tp. We find 



7 MS 



1 + ^C F a 
3e 



+ 
+ 



4 (AT F N f + 8C F - 11C A ) ^ + ^('UCa-UCf -16T F N f )^ 
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C F a 2 



— (l21C| - 132CaCf - 88C A T F N f + 32Cf + 48C F T F N f + lQTpNfj 



+ — [718C A C F - 823Cj + 5UC A T F N f - 168C F 
- U0C F T F N f - 64T|iY/) 



*The symmetrization and traceless symbol, S, is understood in all our subscripts. 
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+ — ( (648C(3) + 2615) C\ 
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(1944C(3) + 1066) C A C F 



(1296C(3) + 782) C A T F N f + (1296C(3) - 70) Cf 



+ (1296C(3) - 853) C F T F N f - IttTpNf) - 



C F a 3 + 0(a 4 ) (2.5) 



which is clearly gauge independent and agrees exactly with the result of 7 . Throughout we will 
present each renormalization constant as an aid to interested readers who wish to perform the 
check that the double and triple poles in e can be derived from the simple poles of the previous 
loop order which follows from the fact that MS is a mass independent renormalization scheme. 
Further, in Q2.5|) C( n ) is the Riemann zeta function and we have set a = g 2 /(167r 2 ), T a T a = C F I, 
jacdjbcd _ (j A §ab anc j ^ ^rparpb^ _ j^ F §ab Consequently, having checked that our programming 
reproduces known results it is elementary to change the renormalization scheme and deduce 
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1 + 



+ 



+ 



C F a 



nc A ) 



+ 



(47C A + (48 + 18a) C F - 16T F N f ) 



(4T F N f + 8C F (2 2 „ 

— ((81a 3 + 486a 2 - 324((3)a + 2439a - 3564C(3) + 12808) C A 

+ (l62a 2 + 630a + 2592C(3) - 456) C F - (720a + 5336) T F N f 



— (l21C| - 132C A C F - 88C A T F N f + 32C F + 48C F T F N f + l&TpNfj -j 



+ 



+ 



\28T F Nf 



— ((413 - 297a)C A C F - 16A6C\ + 1088C A T F N f + (216a + 408)C| 

+ (108a + 92)C F T F N f 

— ((243a 3 + 1458a 2 - 972C(3)a + 9855a - 18468C(3) + 42902) C A C F 

+ (2592C(3) + 10460) C\ - (5184C(3) + 3182) C A T F N f 
+ (486a 2 + 1134a - 4252 + 12960C(3)) C\ 



U8T F Nf 



1 



+ (5184C(3) - 22396 - 3024a) C F T F N f ) - 



+ 



f f 8748a 5 + 83106a 4 - 17496C(3)a 3 + 493776a 3 + 1659204a 2 

69984 VI sw 

- 96228C(3)a 2 - 2523312C(3)a + 174960C(5)a + 7952229a 

- 11944044C(3) + 746496C(4) + 524880C(5) + 38226589) C\ 
+ (l7496a 4 + 95256a 3 + 151632C(3)a 2 + 416016a 2 

+ 346032C(3)a + 466560C(5)a + 921564a - 4914432C(3) 

- 2239488C(4) + 8864640C(5) + 3993332) C A C F 

- (77760a 3 + 466560a 2 - 124416C(3)a + 4091040a 

- 369792C(3) + 1492992C(4) + 24752896) C A T F N f 

- (155520a 2 - 622080C(3)a + 1937088a - 3234816C(3) 

- 1492992C(4) + 9980032) C F T F N f 

+ (345600a + 221184C(3) + 3391744) T F Nf 

+ (l7496a 3 + 289656a 2 - 373248C(3)a 2 - 715392C(3)a 
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+ 879336a + 10737792C(3) + 1492992C(4) 
- 9331200C(5) -3848760)^)1 C F a 3 + 0(a 4 ) 



(2.6) 

It is worth commenting on the status of the variables in the renormalization constants in both 
schemes. We adopt the convention that the scheme is indicated on the renormalization constant 
itself and therefore the coupling constant and covariant gauge parameter are the variables of 
the same scheme. However, the three loop renormalization of QCD in the RI' scheme has been 
given in [32] and the relationship between the parameters has been determined, [32], as 



and 



l Rl' 



MS 



+ 



Ms) 



(2.7) 



a 



RI' 



1 + 



+ 



.0^2 _ 1 

MS MS 



97 C A + S0T F N f 



36 



+ ( ( ISaijg - ISofa + 190a^ - 576C(3)a^ + 463a^ + 864C(3) - 7143) C| 



MS 



MS 



MS 



'MS 



+ 



+ 



320a 



MS 



3200^ + 2304C(3) + 4248 J C A T F N f 
4608C(3) + 5280 J C F T F N f 



(( 



486a Hs + 1944q HS 



2 

288 



+ 4212C(3)a| I g 



5670C(5)a| 



18792a- 



+ 48276C(3)a| I g - 6480C(5)a|^ - 75951a 
+ 2916C(4)a^ + 124740C(5)a^ + 92505a|^ 



MS 



MS 



+ 



MS 

52164C(3)a 

m , ^.^^^ m , - 1303668C(3)a M g 

+ 11664C(4)a M g + 447120C(5)a M g + 3548070^ + 2007504C(3) 
+ 8748C(4) + 1138050C(5) - 10221367) C\ 

(l2960a|^ - 8640a|jg - 129600C(3)a|jg - 147288af^ 
+ 698112C(3)a M g - 3123360^ + 1505088C(3) - 279936C(4) 
- 1658880C(5) + 9236488) C 2 A T F N } 



MS 



+ 



+ 



+ 



(248832C(3)afjg - 285120a|^ + 248832C(3)a M g 



2851200^ 



( 



5156352C(3) + 373248((4) - 2488320C(5) + 9293664) C A C F T F N f 



38400a| 



MS 

884736C(3) 
3068928C(3) + 4976640C(5) 



221184C(3)a M g + 55296a^ 



MS 



1343872) C A T F N? 
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988416) C F T F Nf 



+ (2101248C(3) - 2842368) C F T F Nf 



a- 



MS 



31104 



a MS + 



°Ms) 



(2i 



where the scheme of the variable is indicated as a subscript and we note also that the RI' and 
MS /3-functions are formally equivalent at four loops, |32| . 

From these renormalization constants we can deduce the anomalous dimensions from 



and 



/3(a) 



da 



da 



da 



(2.9) 



(2.10) 
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We find, in our conventions 

„MS 



Mb / \ 



-C F a + 

3 27 



+ 



8C F 
~243 



[47C A - UC F - WTpNf] a 2 
(648C(3) + 2615) C\ - (1944C(3) + 1066) C A C F 



(1296C(3) + 782) C A T F N f + (1296C(3) - 70) Cf 



+ (1296C(3) - 853) C F T F N f - ll2T F Nf\ a 3 + 0(a 4 ) (2.11) 



which agrees with 7 and 



8„ , C F 
3 54 



Cf_ 
1944 



(27a 2 + 81a + 1434) C A - 22AC F - 504T F N f a 2 

(486a 4 + 4131a 3 - 972((3)a 2 + 19899a 2 - 10044C(3)a 

+ 64098a - 115344C(3) + 806800) C 2 A 

- (l458a 2 + 15066a + 10368C(3) + 172856) C A C F 

- (4320a 2 - 2592C(3)a + 22248a + 25920C(3) + 534400) C A T F N f 
+ (3888a + 41472C(3) - 43328) C F T F N f 

+ (82944C (3) - 4480) C 2 F + 782087^/1 a 3 + 0(a 4 ) (2.12) 

in four dimensions. We adopt the same convention for the renormalization group functions as 
for the renormalization constants in that the scheme the variables correspond to is indicated on 
the left hand side. Restricting to SU(3) and the Landau gauge we have 



SU(3) 
a=0 



32 4 

~9~ a " 243 [378Nf " 6005] ° 2 



+ 



6561 



10998JV/ - 6318C(3)iV / - 467148JV/ 
- 524313C(3) + 3691019] a 3 + <3(a 4 ) 



(2.13) 



As an aid to the lattice matching procedure we note that the finite parts of the various compo- 
nents in the two schemes after renormalization are 



y (l) MS finite/ x 



1 - -C F a 



'9487 



1 



a a 

V 324 3 8 



8C(3) + 2C(3)a ) C F C A 



2101 
162 



T F NfC F + [-a -a 



635 452579 976 4 

a rrr^ — C(3) — -C(3)a 



216 



2187 



81 



9 



8195 
648 
64 
T 



+ 8C(3)J C 2 F I a 2 
C(4)) T F N f C F C A 



85 61322 272 w . 

—a + C 3 

9 2187 9 w 



y C(3)a + ^C(4) ) T F N f C 2 



+ 
+ 



/ 63602 256 
V 2187 + "81 
/ 21026833 953 



+ ^C(3) ) T 2 N 2 C F 



75 



V 69984 



1728 a+ 64 a + 48 a 



/135 5 „ . . 
{T + 4° I C(5) 



+ 
+ 



1 



-a 3 + 2a 2 + 



167 
36" 



3230 , 

S-r-^ 2 )« 4 ») c ^ 



7 



+ 



+ 



-a 



41 

72 



a 



979 
~27 



a 



605431 
4374 



/440 40 \ A . r , 



- - f a " \o? + a 3 ) C(3) - 38C(4)) C 2 F C A 

Z7964 164 16 2 2 3 \ 64 400 

a 3 + 0(a 4 ) (2.14) 



2935 821 79 , 

1 a H a 

5832 36 36 



y (2) MS finite, x 



and 



y (2) RI' finite,. 



+ 2a C F a 



284 29 , 
— + T a + 2a^ 

20. 



6C(3) - 2C(3)q C f C a 



35 



74 



-T F N f C F + [a' - —a - — + 16C(3) C 



/184 > . . 1 3 6 i . 
(-§-C(3) + -g-C(3)a 



9 

2267 



+ 
+ 

+ 

+ 



( 



1399 4913 
a + 



54 
1088 



-a 



87980 
243 



T F NjC F C A 



C(3) ) iwy^ + ^n*$c F 



V 81 27 9 * \ * j 1 j r 243 

/ 2762093 28855 431 , 37 , /445 \ „, , 

(^888- + + "16^ + ¥ " + I" + 5a j C(5) 

'23 2 H65 1617\ \ 2 



a' + 



18 



o 55 2 4790 
+ —a - 



4 
5. 
36 
'7450 



-a + 



81 

/7450 131 
+ (^ + ^ 
/1760 W . N /4232 



a — 



4 

198059 
324 



C(3) C F C A 



(^-f a )c(5) 



a 



5a 2 ) C(3)) C 2 F C A 



V 3 



-C(5) 



V 9 



+ 8a K(3) 



3281 79 2 1 3 \ „i 

+ -32T a+ i8 a ~r ]cf 



29483 
~ 972 

a 3 + 0(a 4 ) 



(2.15) 



- + 2a \ C F a 



(284 271 1 \ 

-g- + ^g-a + 3a 2 + -a 3 - 6C(3) - 2C(3)aj C F C A 



'20 



+ 



40 \ 
— a J TWV/Cf 

604 



+ (« 2 + 3a - ^ + 16C(3) J C\ 



a" 



/184 



■C(3) + f C(3)a - ^a 3 
40 , 9470 87980 \ 

m 36536 1088 

—a 2 + 52a - + ■ 



-a z 



+ 



/9680 800 N 



/yoeu ouu \ 2 2r , 
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+ 



+ 



+ 



/ 2762093 289991 



V 3888 



+ 



569 



-a H of + 

1296 12 



127 



19 



1 



11 



—or H a + 

2 « t- 4 « 



685 
~9~ 



1617 



-Q 



9 a + + 4° 



C(3) 



/ 445 
I 3 



+ 

1 4 5 3 

r + 2° + 



+ 5a C(5) C F Ci 



107 



+ 



/7216 



V 9 
1760 



V 3 



C(5) 



9 

7a — a 

Z3688 
\9~ 



a' + 

2 



1429 



469555 



-a 



36 972 
/ 1480 40 



C(3) 

8a] C(3) 



35 



31 



1 



— a H a H — a" 

6 6 2 



C 3 



56287 
486 

: + 0(a 4 ) 



a C(5) C\C A 



(2.16) 



These imply 



y (i) MS finite, v 



SU(3),a = 



124 

451293899 
157464 

( 



( 68993 160 



V 729 9 
1105768 



C(3) 



2101 
243 



N, 



2187 



C(3) + ^ ((4 ) + ^C(5) 



/ 63602 



324 
8636998 



81 



19683 



V 6561 



+ *? 



+ 0(a 4 ) 



(2.17) 



and 



y (2) MS finite, x 
ip^D v tp W> 



SU(3),a = 



9° 



Z2224 40 w , 40 , T 



136281133 
26244 
/1232 



S7M4I «rao 



243 



+ 0(a 



+ 



V 81 



C(3) 



15184 
27 



81 

»r 9680 »- 

iV/ + -^-iv; 



729 



(2.18) 



3 Second moment of transversity operator. 

The renormalization of the transversity operator ipa^ 'D p ip is similar to that for ■07 /J D h 'ip aside 
from the decomposition of the corresponding Green's function. The specific operator we are 
interested in is 

(d-1) 

+ 7^31) (v^^Dx^ + v^v uX D x 4,) (3.1) 
which is symmetric with respect to the indices v and p and satisfies the traceless conditions, 

Sipa^Dytl) = S^o^Dvi) = . (3.2) 
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The Lorentz decomposition of the particular Green's function we will renormalize is 

~,pyp , 



(p) = m P ) [S^D^](0) ^(-p)> 



(p) (a^pO + a^p v - o^pPpx + ^a^px^j 



+ E 



(2) 

i\)<jP v DP i\) 



(p) ( rf v aP x px + if p a vX px + ^ 1 ^ + 2 ^ cj^Vp p Pa 



- (d+l^^vV) 



(3.3) 



where the components are determined through 



ipaP v DP'4> 



(p) 



8(d- l)(d-2 



-tr 



(d + 2) A 

0>Pp + 0>p£^ 2 a ^PvPpP 



8(d- l){d-2)p 2 



tr 



{o v xPpP P p X + v P \PtiPuP X 



+ da ii xPvP P p X ~ r} up a^ x p X p 2 ) G^ v {p) 



v( 2 ) - 



(d 2 - d + 2) 



tr 



(j]pv°p\P X + ri w a v xp X - (d + l)v^pXP X 
(d-l)(d + 2) 



+ 



VpXPuPpP )G^ vDPip (p) 



+ 



,(3) 



(P) = - 



8(d + l)(d- 2)p : 



tr 



■tr 



{o v xPp,PpP X + (yp\Pp,PvP y 



+ 



tr 



8(d-l)(d-2) 
1 

8(d + l)(d-2) 

(d + 2) 
8(d + l)(d-2)p 2 



+ da^xPvPpP -Vupcr^p p J G^ DP ^(p) 
(d + 2) \ 

0>Pp + CTppPz. ^ CTpXPuPpP 



+ VupVpxP ) G^ DP ^(p) 

(vpis (J p\P X + Vpp°v\P X - (d + l)ri up anxP 
(d- l)(d + 2) 



+ 



tr 



(<Tu\PpPpP X + cr pX PpPuP X 

+ da flX PuPpP X ~ Vvp(Tp\P X p 2 ) G ^ VDP ^{p) 



(3.4) 



Using the same programmes which determined the renormalization of ip^D v ip but instead 
including the Feynman rule for the transversity operator, we find after renormalizing the Green's 
function in the MS scheme that 

7 MS - 3 

Zj tpaP v DPip 



I -'- -C F a 
e 



+ 
+ 



(4T F N f + 9C F -llC A )^ + (35C A -9C F -l2T F N f )±- 



C F a 2 



(242C\ - 297C A C F - 176C A T F N f + 81Cf + 108C F T F N f + 32T|iV / 2 ) 



18e 3 
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+ (lU3C A C F - 1178C1 + 784C A T F N f 
- 243C 2 F - 252C F T F N f - 96T F N^ 
+ (l2553Cl - 7A79C A C F - (5184C(3) + 4168) C A T F N f + 1782C| 



+ (5184C(3) - 3240) C F T F N f - S68T F Nf 



324e 



C F a 3 + 0(a 4 



and with the RI' scheme definition 



lim 



(P) 



= 1 



(3.5) 



(3.6) 



we have 



7 RI' 



1 + 

+ 



+ 



- e +\(3a + 7) 



C F a 



1 



(4T F Nf + 9C F - IIC A ) 7^2 + ( 35C ^ + (33 + 18a)C F - l2T F N f ) 

+ — ((9a 3 + 54a 2 - 36C(3)a + 271a - 324((3) + 1189) C A 

+ (36a 2 + 108a + 288C(3) - 165) C F - (80a + 468) TpiV/)] a 2 

(242C1 - 297C A C F - l76C A T F N f + 81C F + W8C F T F N f + 32T|iV / 2 ) ^ 

+ ((450 - 297a)C j4 C F - 1178(7^ + 784C A T F N f + (243a + 324)C|. 

+ 108aC F T F N f - 96T 2 iV/) J-^ 

((729a 3 + 4374a 2 - 2916C(3)a + 30456a - 26244((3) + 101196) C A C F 
+ 25106C1 - (10368C(3) + 8336) C A T F N f 
+ (2916a 2 + 6561a - 14904 + 23328C(3)) C 2 F - 736T F Nj? 

+ (10368C(3) - 51192 - 9396a) C F T F N f ) 



+ 



648e 



+ 



5184 



— ((972a 5 + 9234a 4 - 1944C(3)a 3 + 54864a 3 - 10692C(3)a 2 

+ 184356a 2 - 282312C(3)a + 19440C(5)a + 874185a 

- 1644228C(3) + 447120C(5) + 3930425) C\ 
+ (3888a 4 + 24300a 3 + 3888C(3)a 2 + 116748a 2 

- 34992C(3)a + 51840C(5)a + 375732a + 1260576C(3) 

- 466560C(5) - 511188) C A C F 

- (8640a 3 + 51840a 2 - 13824C(3)a + 454560a 

- 124992C(3) + 124416C(4) + 2438256) C A T F N f 
+ (38400a + 18432C(3) + 323200) T F Nf 

- (34560a 2 - 82944C(3)a + 295920a - 82944((3) 

- 124416C(4) + 620592) C F T F N f 
+ (7776a 3 - 31104C(3)a 2 + 54432a 2 + 32076a + 51840C(3) 
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+ 207360C(5) -266868) Cf) C F a 3 + (9(a 4 ) 



(3.7) 



Consequently, the renormalization group function 



^ S (a) 



3C F a + ^f[35C A -9C F -l2T F N f ] 



+ 



Cf_ 
108 



12553C1 - 7479C A C F - (5184C(3) + 4168) C A T F N f 
+ 1782C| + (5184C(3) - 3240) C F T F N f - 368T F Nj> 



+ 0(a 4 ) 



(3.8) 



has emerged as being independent of a as expected. The two loop part of (j3.8|) agrees with the 
known results, jHSl EM IHB Ell SU , when they are restricted to the second moment. A second 
check on the result is that the double and triple poles in e at three loops in the renormalization 
constant have correctly emerged. Also, in four dimensions 



_RI / \ 



3C F a + ^ 



+ 



Cf_ 
432 



(9a 2 + 27a + 364) C A - 54C F - l28T F N f ] a 2 

(l62a 4 + 1377a 3 - 324((3)a 2 + 6633a 2 - 3348C(3)a 

+ 21366a - 42768C(3) + 224296) C\ 
+ (l62a 3 + 324a 2 - 1674a + 38016C(3) - 71100) C A C F 

- (l440a 2 - 864C(3)a + 7416a + 5184((3) + 145600) C A T F N f 
+ (432a + 6912C(3) - 4032) C F T F N f 

+ 7128C£ + 20992T|JV^1 a 3 + 0(a 4 ) . (3.9) 
Thus, for SU(3) and the Landau gauge we have 



-y^ s (a) 



SU(3) 
a=0 



4o - 2 [2N f - 31] a 2 



81 



and 



Til' (n \ 



SU(3) 
a=0 



4a 



92iV/ + 4320C(3)A / + 8412A^ - 86229 a 3 + 0(a 4 ) (3.10) 



[16JV) - 255] a 2 



+ 81 
+ 0(a 4 



656JV? - (396C(3) + 27636) N f - 29106C(3) + 218367 



(3.11) 



which are the main results of this article. To assist with the lattice matching procedure we record 
that the finite parts of the components in both the MS and RI' schemes after renormalization 
are 



y (i) MS finite/ n 



1 



7 «w> 
— + — ) C F a 



2 

/943 3 3 2 

1 — a H — a 

V 24 4 8 

37 5 3 5 
T + 2 a "4 a 



^C(3) + ^C(3)a ) C F C A - 16T F N f C F 



12C(3) C F 



12 



and 



/478 47 \ 

\27 + 1° " 8C(3) " + 24C(4) J TpN f C F 



+ 
+ 



+ 



+ 



(2) MS finite, x 



y (3) MS finite, . 



1160 



V 27 



+ f C(3) ) WCfc 
61 10 



y + y«)C(3)-24C(4) 



104843 68 



324 

/ 2656369 18379 253 2 113 , 

1 a H a H a" 

V 5184 576 32 96 



) T F N f C F C A 



( 26839 131 



1 a 

144 12 



16 



-a 



^a 3 ) C(3) 



69 
16 



3 
16 



« 2 ) C(4) 



/1699 131 
- 6C(4) 



-a 



4 4 
158969 2657 



a 3 C(3) - (70 + 10a) C(5) 



432 



48 



a - -a 2 + -a 3 C 2 F C A 



40C(5) 

521 59 
12~ + ~2 
1 



74 + 24a + 6a 2 - -a 3 j C(3) 
3 



33 



+ — + —a + -3- a - o a C'f 
8 8 / 

(1) MS' finite, ' 
3 ftj^DPtp W 



a 3 + 0(a 4 



r (3) RI' finite, v 



to three loops. Hence, 



MS finite, x 



5C/(3),a = 



14 

— — a — 
3 

"1852993 


Z2237 
V 18 
97391 


432 


108 



+ 0(a 4 ) 
= 0(a 4 ) 



, 79 A/ „ 7060 
C(3) + T C(4) + ^- 



(3.12) 
(3.13) 

(3.14) 



C(5) 



+ 



/122 

1160 
81 



C(3) - ?C(4) 



306881 
486 



+ §c(3))iv; 



a 3 + 0(a 4 



(3.15) 



4 Conversion functions. 



As in [32 we have checked our expressions by also constructing the conversion functions Co (a, a) 
for each of the operators explicitly. These functions allow one to move from one scheme to 
another and are defined by, [l2*] . 

z m' 

C (a,a) = J=. (4.1) 
The renormalization group functions are then related by 

TO 1 ' Hi') = lo S ( a Ms) " P ( a Ms) ln ( a MS' a MS 

oa MS 
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where we have indicated the scheme of the variables explicitly. Being careful to first convert the 
variables to the same reference scheme, which we will take to be the MS scheme, we have from 
the various renormalization constants, 

C F a 

C^ D „^(a,a) = 1 + (9a + 31 



+ 



+ 



+ (81a 2 + 315a + 1296C(3) - 228) C F - 2668T F N f 
(l61838a 3 - 201204C(3)a 2 + 942597a 2 - 2124792C(3) 



9 

(l62a 2 - 162C(3)a + 783a - 1782C(3) + 6404) C A 

C F a 2 
162 

a 

+ 174960C(5)a + 4796982a - 11944044C(3) + 746496C(4) 
+ 524880C(5) + 38226589) C\ 
+ (26244a 3 + 151632C(3)a 2 + 159408a 2 + 346032C(3)a 
+ 466560C(5)a + 554904a - 4914432C(3) 

- 2239488C(4) + 8864640C(5) + 3993332) C A C F 
+ (528768C(3)a - 1469016a + 369792C(3) 

- 1492992C(4) - 24752896) C A T F N f 

+ (l7496a 3 - 373248C(3)a 2 + 289656a 2 - 715392C(3)a 
+ 879336a + 10737792C(3) + 1492992C(4) 

- 9331200C(5) - 3848760) C 2 F 

- (497664C(3)a + 351648a - 3234816C(3) 

- 1492992C(4) + 9980032) C F T F N f 

+ (221184C(3) +3391744) T 2 iV/j + 0(a 4 ) (4.3) 



and 



C^dp^(«>«) = 1 + (3a + 7) C '' 



+ 



+ 



+ 



(36a 2 + 108a + 288C(3) - 165) C F - A86T F N f 



2 

(36a 2 - 36C(3)a + 174a - 324((3) + 1189) C A 

C F a 2 
24 

(l7982a 3 - 22356C(3)a 2 + 104733a 2 - 238032C(3)a + 19440C(5)a 
+ 523602a - 1644228C(3) + 447120C(5) + 3930425) C 2 A 

+ (l0692a 3 + 3888C(3)a 2 + 63180a 2 - 34992C(3)a 
+ 51840C(5)a + 312876a + 1260576C(3) 

- 466560C(5) - 511188) C A C F 

+ (58752C(3)a - 163224a + 124992C(3) 

- 124416C(4) - 2438256) C A T F N f 

+ (7776a 3 - 31104C(3)a 2 + 54432a 2 + 32076a 
+ 51840C(3) + 207360C(5) - 266868) C\ 
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- (41472C(3)a + 101520a - 82944C(3) 
- 124416C(4) + 620592) C F T F N f 

+ (18432C(3) + 323200) T F Nf\ + 0(a 4 ) . (4.4) 

With these expressions given in terms of the MS coupling constant and covariant gauge param- 
eter, we have checked that the same RI' scheme anomalous dimensions emerge for each operator 
and for non-zero a. 



5 Large Nf critical exponent for transversity operator. 

In order to provide another independent check on our transversity results for MS we have ex- 
tended the leading order large Nf calculation of the critical exponent associated with the flavour 
non-singlet twist-2 operator anomalous dimension as a function of d = 2\x and arbitrary moment 
n, [331 j to the transversity case. In other words we consider the operator ipcr^ 1 D 1 " 2 . . . D^-^ip 
which is symmetric in the indices {v{\ and satisfies a more general traceless condition to that for 
n = 2, [""7]. As the large Nf critical point technique has been widely documented in [HI Hoi H5] 
we refer interested readers to these sources and quote the main result of our computations. Fol- 
lowing the procedure of [13], we found the <i-dimensional critical exponent for the transversity 
operator is 



V 



(n) 



2(/x - l)C F rft 
(2fi-l)(n-2)T F N f 



O-i) 



Ufa - 2f 



(ji + n- l)(/i + n-2) 



+ 2fi + n - 1) - ^(/i)] 



+ O 



(5.1) 



where 



(2 M -l)(/x-2)r(2/x) 



4r2(^)r(^ + i)r(2- / ,) (5 ' 2) 

and ip{x) = d\n.T(x) / 'dx and T(x) is the Euler gamma function. If one expands (|5.1[) in powers 
of e, [i = 2 — e, then there is a one-to-one correspondence with the leading large Nf coefficients 
of the associated anomalous dimension to all orders in the (perturbative) coupling constant. 
Moreover, since we have a result for arbitrary n we can deduce an analytic expression for the 

leading large Nf part of the three loop term of l^f^ DP ^{ a ) an d check that it agrees with our 
result for n = 2. In particular if we formally write 



7 W MS ( a) 



c F 



oo r— 1 



ha + (b 21 T F N f + b 20 ) a 2 + J2J2 b rj T F N j 



r=3 j=0 



(5.3) 



where the coefficients, {bij}, are functions of n and bi, 620 and 621 are known, [3*5113*9*11101 13TI 11*1] . 
then we find the leading order three loop coefficient is 



4 
27 



48S 3 (n) - 80S 2 (n) - 165i(n) + 



3(17n 2 + 17n 



n(n + 1) 



(5.4) 



where Sj(n) = Ei=i l/« • Setting n = 2 we find exact agreement with the corresponding term of 
(|3.8|) . Moreover, we can expand the critical exponent to four loops and deduce, using the four 
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loop /3-function, [ifi] . 

6 43 = _ [965 4 (n) - 160S 3 (n) - 325 2 (n) + 192C(3)5i(n) - S2S 1 (n) 
81 |_ 

, AA „ n . (131n 4 + 262n 3 + 211n 2 -16n-48)l 

- 144C3 +" 2( , U2 • 5.5 

n z {n + l) z 

6 Discussion. 

To conclude we have derived expressions for the three loop anomalous dimension of the second 
moment of the transversity operator in both the MS and RI' schemes. The results for the latter 
will be important in the extraction of lattice estimates of the associated matrix elements. More- 
over, we have provided the same information for the usual twist- 2 flavour non-singlet operator. 
Although we have had to consider various projections of the Green's function with the operator 
inserted in a quark two-point function it ought to be possible to extend the present calculations 
to derive the three loop anomalous dimensions for higher moments of the transversity operator 
in the MS scheme. In this case since one is only interested in the divergent part and not the 
finite part of the Green's function, it would not require the same use of the projections consid- 
ered here. This is important since the algebraic manipulations are slowed for a large number of 
projections and the inclusion of more covariant derivatives in the operator itself. However, since 
those anomalous dimensions are independent of the gauge parameter, the algebraic computations 
would be speeded by considering the Feynman gauge. 
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